#### Generate posterior predictive distribution draws
library(uds)

# For paper

load("post.uds.RData")
ppd <- pred.y(post.uds)
save(ppd, file="ppd.RData")
rm(list=ls()) # free up some RAM

# Robustness stuff

#load("post.decile.RData")
#ppd <- pred.y(post.decile)
#save(ppd, file="ppd.decile.RData")
#rm(list=ls()) # free up some RAM

#load("post.vintile.RData")
#ppd <- pred.y(post.vintile)
#save(ppd, file="ppd.vintile.RData")
#rm(list=ls()) # free up some RAM
